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Abstract 

We present a formal wave theory for the calculation of the spectrum and the eigenmodes for 
a certam class of ray-chaotic optical cavities introduced by A. Aiello, M. P. van Exter, and J. P. 
Woerdman |quant-ph/0307119| . 
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In a previous paper Q], we presented a theoretical model for a composite optical cav- 
ity made of standard laser mirrors; the cavity consists of a suitable combination of stable 
and unstable cavities as shown in Fig. 1. By using numerical simulation we were able to 
demonstrate that such a cavity displays classical (ray) chaos, which may be either soft or 
hard, depending on the cavity configuration. In this paper we want to go a step further 
by addressing the behavior of the chaotic cavity in a wave regime (or, loosely speaking, in 
a "quantum" regime Q]). More precisely, in this paper we present a formal theory for two 
coupled unstable cavities. We show that it is possible to introduce an unitary coupling which 
accounts both for direct transmission and diffraction (which occurs from the edges of the 
convex mirrors in our cavity) by using a suitable scattering operator (see Eqs. (0121) below). 

A standard two-mirror stable resonator is a geometrically open system but because of its 
stability it is closed both from ray B] and wave point of view. In other words, a typical 
gaussian-beam-like mode in such a resonator is confined both longitudinally (that is along 
the axis of the resonator) and transversally (that is along the two directions orthogonal to 
the axis) by the focussing action of the two mirrors. Because of this confinement a stable 
resonator has a discrete spectrum; in paraxial approximation this spectrum can be classified 
in a "longitudinal" part which depends only on the length of the cavity and in a "transver- 
sal" one which depends also from the radii of curvature of the two mirrors. Here we are 
interested mainly in the transversal part. 

Efficient methods to calculate the spectrum and the eigenmodes of hard-edged unstable 
cavities were developed in the last 30 years; particularly notable is the asymptotic theory 
created by Horwitz and Southwell p. However, in spite of this long hystory, surpris- 
ing properties of these eigenmodes were discovered recently 0, Q, ^- For instance, the 
Horowitz-Southwell theory has been exploited and slightly modified by Berry et. al. to 



investigate both the fractal nature of the cavity eigenmodes j9| and the occurrence of the 

n 

Petermann excess- noise factor 10']. In this paper we apply Berry's theory to our composite 
cavity, thus generalizing some of the results presented in jsl. From a mathematical point 
of view, the main difference between the theory for a conventional unstable cavity and our 
composite system, is that in the former case the operator which accounts for the modes 
propagation inside the unstable cavity is not unitary because of the losses from the edges 
of the smallest mirror. As we shall show later, in our case the two round-trip operators de- 
scribing the mode propagation in the two half cavities shown in Fig. 1 remain non-unitary 
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but the operator describing the motion in the overall cavity is unitary because the whole 
cavity is stable {L < 2R). 

In this paper we restrict our attention to two-dimensional cavities with one-dimensional 



mirrors {strip resonators). Following Berry li| it is convenient to introduce from the begin- 



ning a "quantum-like" vector-space notation writing the modes of the field as kets in a linear 
space defined by the propagation operator K whose coordinate representation is given by the 
Huygens' integral in the Fresnel approximation Within this formalism, the transversal 
mode profile u{y) calculated in an arbitrary plane z = const, can be considered as the coor- 
dinate representation of a field state \u) depending on the longitudinal coordinate z which 
is considered as a parameter (exactly as the time in the Scrodinger equation): 

{y\u) = u{y). (1) 

In order to describe the dynamics of each sub-cavity and the coupling between them, we 
introduce a set of four fields Ui,U2 and fi,f2 defined in the reference plane z = following 
the scheme illustrated in Fig. 2. Then the propagation in the left and right side of the whole 
cavity can be described by introducing the operators Kl and Kr respectively: 

= e * A Kl\v2}, 
\U2) = e * A Kr\vi}. 



At this point the two sub-cavities are still uncoupled. In Eq. Q Kl = K{li), Kr = K{h), 
where /i and I3 are the lengths of the left and right cavity respectively and the coordinate 



representation of the paraxial propagator is 



{y\m\y') = V^^^p [-^ (^^" " '^yy' + ^^')] • 

The three coefficients A, D, B are the corresponding elements of the following ABCD matrix: 



M(0 = ^ ^ ' f \. (4) 




where A = D. 

In order to describe the coupling between the two half cavities we introduce the four 
scattering operators Sij = 1,2) 

\vi) = Sii\ui) + 5*121^2), 

|^^2) = 5'2i|Mi) + S22\U2), 



3 



where the diagonal operators Sa describe the transmission of the field above the central 
mirror {\y\ > a) while the off-diagonal operators Sij {i j) describe the reflection on the 
central mirror {\y\ < a). We require that the coupling between the two half cavities is 
unitary by imposing: 

+ (M2IM2) = {Vl\Vl) + {V2\V2), (6) 

from which it follows that: 
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Y,SlS,k = iS^k, {^,J,k = 1,2), (7) 

where Sik is the Kroneker tensor. Since the bi-convex optical element in the center of our 
cavity (see Fig. 1) is invariant with respect to the symmetry z —z, we can assume that 
the coupling is the same going from left to right and viceversa, and put: 

5*11 = 5*22 = T, S12 = S21 = R, (8) 

from which it follows that the unitarity conditions Eq. ((7j) become: 

f^f + B)R = 1, 

.... (9) 
T^R + R^T = 0. 

Before investigating the consequences of these relations we collect the four fields Ui,U2 and 
f 1 , V2 in doublets 

{ui,U2}^ I '""'^ I , {v^,V2}^ I '""'^ I , (10) 

\ 1^2) J \ \V2) J 

which represent the incoming and outgoing fields in the plane z = respectively. Alterna- 
tively is possible to relate the fields in the left side of the cavity {ui,V2} with the fields on 
the right side {^1,^2} by introducing a set of four transmission operators that are related 
in a simple way use the scattering 

formalism. Now we can rearrange the previous Eqs.(j2]ini) 






Inlo . 




and 





(11) 



(12) 



respectively. Inserting Eq. ()12|) in Eq. we obtain, after a few straightforward algebraic 
manipulation, and assuming the simpler case li = I3 = I ^ Kr = Kl = K, the eigenvalue 
equation for the modes of the cavity: 

M-7i TK \ ( \v,)\ ^ ^ 

TK RK--fl J y \v2) J 

where we defined the eigenvalue 7 as: 7 = exp(i^). By inspecting Eq. (fT^ we can easily 
recognize that the product RK = Krt is the well known round-trip propagator for a 
single sub-cavity. Moreover we notice that when T = we get two independent eigenvalue 
equations for the two unstable sub-cavities; in this case Krt is not longer unitary and 
I7I < 1. With Eq. (fT!^ we have achieved the goal of this paper. This equation can either 
be solved numerically by diagonalizing the matrix in Eq. ()13|) or by applying asymptotic 
methods jsl. 

In order to write Eq. ()13p in coordinate representation is necessary to write down the 
explicit form for the transmission T and the reflection R operators. To this end we first 
notice that the paraxial propagator which accounts for the reflection by a convex mirror has 
the following coordinate representation: 



/ 2TTi \ 

{ylfW) =expi-—y^j5{y-y'), (14) 



where r is the radius of the convex mirror Since the reflection operator is a mathemat- 
ical representation of the bi-convex mirror whose transverse dimension is 2a, its coordinate 
representation must be limited to the region \y\ < a. Analogously it is easy to understand 
that the transmission operator can only exists in the region \y\ > a. These physical consid- 
erations make it natural to try the following expressions for the transmission and reflection 
operators: 

{y\f\y') = 5{y-y')Q{\y\-a), 

{y\R\y') = 6iy-y')Qia-\y\)exp{-'^y'). 

It is easy to check, by straightforward calculation, that choosing this form for the R and T 
operators, Eqs. Q are automatically satisfied because of the following properties of the 9 
functions: 

Q(\y\-a) + eia-\y\) = 1, ^^^^ 
e{\y\-a)eia-\y\) = 0. 
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In conclusion, we have derived the equations for a pair of coupled unstable cavities. We 
obtained an eigenvalue equation ()13|) which can be solved in straightforward way to get the 
spectrum and the eigenmodes of the whole cavity. The theory in the present form involves 
some not well defined quantities (as products of distribution functions) which are justified 
only on a physical basis. 

This project is part of the program of FOM and is also supported by the EU under the 
IST-ATESIT contract. 



APPENDIX 



In this appendix we give some details about practical calculations. We start rewriting 
Eq. (0 as 

= {y\RKW) + {y\fk\v2), 

(17) 

= {y\fk\v,) + {y\Rk\v2). 
For simplicity we define RK = p and TK = f and write explicitly Eq. ()17|1 as: 

iviiy) = / dy'p{y,y')vi{y') + / dy'T{y,y')v2{y'), 



(18) 

7V2{y) = I dy'T{y,y')vi{y') + I dy' p{y,y')v2{y'), 



where we have defined p{y,y') = {y\p\y') and T{y,y') = {y\r\y'). For a symmetrical cavity 
with li = /s and we look for a solution such that vi{y) = V2{y), therfore Eqs. (fTHj) reduce to 
a single equation 



iMy) = j dy'[p{y,y') + r{y,y')]vi{y'), 

(19) 

= [e(a-|i/|)e(l/) + e(|y|-a)] [ dy'Kiy,y')vi{y'), 



where we have defined ^(y) = exp Here K{y,y') is the propagator from a round- 

trip inside one unstable sub-cavity without accounting for the reflection on the convex mirror. 
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Instead the product ^{y)K{y, y') = KjiT^y, y') gives us the propagator for a complete round- 
trip. For computational reasons is more convenient to work with RRTiViy') instead of 
K{y,y') therefore, exploiting the fact that = 1 we rewrite Eq. (fTUI) as 



^v,{y) = [0(a - \y\) + e(|y| - a)C{y)] / dy'KnT{y,y>iiy'). (20) 



After scaling all lengths with a, Eq. (p]j)can be written as 



^g{y) = a/- [6(1 -\y\) + Qi\y\ - l)Ciy)] I e-^'^^-y/^'^'' g{x)dx, (21) 



where, following Horwitz Q], we have defined: 



tR 



2l\{l-l/R) ' 

t = ttMF, (22) 

g{y) = e'^FiM-M-^)/2y^^^yy 

The magnification M can be also written in term of m = {A + D)/2, the half of the trace of 
the ABCD matrix, as M = m + V — 1. In practice we have to calculate the asymptotic 
form of the following three integrals: 



oo 

Ji = / e-''^^-y/^'^"g{x)dx, 



1 
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h = j e-^*(^-^/*^) g{x)dx, (23) 



\-''^^-y/^''^"g{x)dx. 



The value x = y/M (with M > 1) for which the phase is stationary can be inside or outside 
the domain of integration depending on the value of y as illustrated in the following table: 
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TABLE I: The real axis (— oo < y < cxd) has been divided in five subsets. For each of them 
the letters Y/N indicate if the stationary point is contained/not contained within the domain of 
integration of the integrals /i,/2 and Is. 
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FIG. 1: Schematic diagram of the cavity model. Two unstable cavities are coupled to form a 
single cavity which is globally stable for L < 2R. The two sub-cavities are unstable for I < R — r 
and stable for R — r < I < R. 




FIG. 2: Logical scheme of the propagation process and of the coupling between the two sub- 
cavities. The dashed line represent the plane z = where the bi-convex mirror is located. Kl, 
Kji are the operators describing the field propagation in the left and right side of the whole cavity 
while T and R describe the coupling between the two sub-cavities. 
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